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Abstract 

We discuss how the electrons in relativistic GRB shocks can reach near-equipartition 
in energy with the protons. We emphasize the non-Fermi origin of such acceleration. 
We argue that the dynamics of the electrons in the foreshock region and at the 
shock front plays an important role. We also demonstrate that PIC simulations can 
directly probe this physics in the regimes relevant to GRBs. 
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Electron acceleration/heating — There is a lore that charged particles 
are accelerated at shocks by the Fermi mechanism. Numerical simulations 
show that although Fermi accceleration may be present, it cannot heat the 
bulk electrons to near-equipartition with the protons, i.e. ee<l. We suggest 
alternative mechanisms that may be at work in relativistic coUisionless shocks. 



Magnetic fields are genera ted at shocks by the Weibel instability (iMedvedev &: Loebl . 



19991 : IChang et al.l . 120071 ). that saturates its linear phase at a relatively low 
magnetic field, near the equipartition with the electrons. At such low fields, 
protons keep streaming in current filaments, whereas the electrons, being much 
lighter than the protons, are quickly isotropized in the random fields and form 
a uniform background. Nonlinear evolution of the filaments leads to further 
ampUfication of the magnetic field up to ~ 0.1 on average (and ~ 1 lo- 
cally in clumps) as one approaches the main shock compression. The average 
Lorentz factor of the electrons is also gradually increasing tow ard the shock 



and eg becomes 30% to 50% around and after the shock jump (IChang et al. 



20071 : Spitkovsky, in prep.). This electron heating we are trying to understand. 



Electrostatic model — The current filaments are formed by the protons moving 
roughly at the speed of light (their Lorentz factor is ~ F). Hence, they are the 
sources of both the magnetic and electrostatic fields. In the absence of strong 
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Fig. 1. Acceleration mechanisms, (a) — the electrostatic model, (b) — the filament 
merger as an example of the induction model. 

electrostatic shielding (simulations seem to show this) E and B fields are re- 
lated as B<E. An electron, moving toward a filament, see Fig.[T^, gains energy 
Ue — eEl ~ eBl, where / is the radial distance the electron travels through 
the strong E field region. We normalize it to the scaracteristic scale in the 
system, c/upp (the typical size of the filaments is few c/upp) as / ~ \{c/ujpp), 
where Upp = (47re^n/mp7p)^/^ is the relativistic proton plasma frequency and 
7p ~ r. The parameter A accounts for the actual geometry of the filaments, 
the electrostatic shielding in plasmas, the effects of the electrons on the cur- 
rent distribution, etc. If most of the electrons reside near the filaments (as 
simulations also seem to show), the electron energy density is estimated as 
Ue = nue — neBXc/ujpp, or in dimensionless units: — A^ee. Although this 
mechanism does not provide the net acceleration (the fields considered are 
potential), the electrons gain energy locally, and can radiate it away if the 
cooling timescale in filaments is short enough, see below. 

Induction model — The current filaments exhibit violent dynamics in which 
the field configuration can change rapidly; these include the filament multiple 
mergers (in the foreshock) and break-up (mostly at the shock jump). Let's con- 
sider a merger as an example, Fig. [lb. Two filaments with the typical magnetic 
field B approaching each other with a velocity, f ~ c, induce the non-potential 
electric field i^ind in-between. The typical size of the region with this field is 
of order the filament size, hence it is (i ~ 6{c/ujpp), where 5 is a dimension- 
less size of the merger region. An electron traversing the region gains energy 
Ue ~ eEi^^d ~ e{v / c)B5c/ujpp. The corresponding electron energy density is 
UeU. Thus, we obtain ee — {v / c)5^/eBi which recovers the previous result once 
A is replaced with {v/c)5. Perhaps not all plasma goes through such regions, 
so this mechnism may explain acceleration of a smaller number of energetic 
electrons from the tail of the distribution. Unlike the electrostatic acceleration, 
the inductive acceleration is "permanent", meaning that the electrons remain 
energetic in the downstream, where the current filaments are essentially gone. 

Relevance to GRBs — The radiative efficiency of a shock is determined 
by how fast the bulk electrons lose their energy via radiation. If the electron 
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synchrotron cooling time is smaller than or comparable to the electron resi- 
dence time in the high field region (near the shock front), this electron will 
radiate away energy comparable to its kinetic energy. The shock will be ra- 
diatively efficient in this regime regardless of the field dynamics in the far 
downstream. The dimensionless cooling time is defined as Tcooi = tcooi^pp = 

( ar-y^B'^ ) (tw^) ' "where tcooi is the synchrotron cooling time, n' is the par- 
ticle density behind the shock measured in the downstream frame, F is the 
shock Lorentz factor. The region of strong magnetic field at the shock front, 
where ~ 0.1 — 0.05 is of the size ds ~ 50c/upp or so. Since the shock 
moves at f = c/3 in the downstream frame, the residence time in the region 
of high field is tres ~ rfs/f ~ (150 — 300)ti;pp^ This estimate, does not account 
for the filling factor of magnetic inhomogeneities, which shortens the effective 
tj-es, and the electron trapping in high- field clumps, which is increasing tj-es- 



If ^cooi^ires, then the electrons lose their energy quickly near the shock jump, 
hence the radiative efficiency is high and such a shock can be seen as a GRB. 
We refer to a shock as the "radiative shock" if Tcooi^300 and as the "weakly 
radiative shock" otherwise (its efficiency depends on the fields in the far down- 
stream, which have not yet been adequately probed in PIC simulations). In 
an extreme case, Tcooi < 1, called the "radiative foreshock" regime, radiative 
cooling will be substantial even before the main shock compression, hence 
cooling may change the entire shock structure. 



Using the standard shock model (iRees &: Meszarod . l2005l ). the comoving den- 
sity behind an internal shock at a radial distance R from the central engine is 
n = ^TiL/ (AirR'^r'^mpC^), where L is the kinetic luminosity, is F of an inter- 
nal shock. The magnetic field and the electron bulk Lorentz factor are fractions 
tB and eg of the post-shock thermal energy density B' = {STTTirrLpC^n'eB)^^'^ , 
7e = (mp/me)Fjee. Parameters and eg are ~ 10% and ~ 50% as follows 
from simulations. Finally, the dimensionless cooling time in baryon-dominated 



internal shocks becomes T, 
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Ri2^B^e ^1 Similarly, we eval- 



uate the cooling time for the external (afterglow) sho cks for the co r istant 



density ISM and the Wind (n oc -R ^) outfiow models (iGranot et al.l . Il999 



Chevalier fc Lil . 120001 ): the details are in (Medvedev & Spitkovsky, in prepa- 
ration). 



The results are shown in Fig. [2l One can see that internal shocks (Fj ~ 4) 
with strong foreshock emission, Tcooi^l, can occur for F<60 and at -R<few x 
10^° — 10^^, well below the photosphere at such low F's. Strongly radiative 
shocks Tcooi ~ 100 — 300 can occur above the baryonic photosphere in a 
relatively narrow, but very natural range of parameters, F ~ 150—400, Fj>2.5, 
Ri ~ lO^^^'^^'^cm and the outfiow kinetic luminosity L ~ lO^^^^erg/s; hence 
they are hkely observable. Since Tcooi oc F~'^, the region of the parameters 
widens greatly with increasing Fj. 
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Fig. 2. (left) — Contours of Tcooi vs. T for the internal shocks for Tcooi = 1, 100, 300. 
Red fihed regions correspond to Tcooi < 1 (radiative foreshock), dark and hght blue 
regions correspond to 1 < Tcooi < 100 and 100 < Tcooi < 300, respectively (radiative 
shock), and the white region corresponds to Tcooi > 300 (weakly radiative shock). 
Here Fj = 4, L52 = 1, = O.IL, t„,_4 = 1. We also mark the radii beyond 
which the internal shocks can form (Ri), the medium is optically thin to Thompson 
scattering (Rph) and the optical depth due to e^-pairs is below unity (i?pair)- {right) 
— Cooling time in afterglows vs. time after the burst, for the ISM and Wind models. 
We use, E = W^"^ erg, = 10 and nisM = 100 cm and a typical z = 2. 

For external shocks, we also see that, except for the very early times, the af- 
terglow emission should be coming from far downstream, not from the main 
shock compression region. However, in the Wind model, the external shock 
can be radiative up to ~ 100 s after the burst, whereas for the ISM model, the 
radiative shock regime can occur only at times earlier than a second after the 
explosion. Since the afterglow usually sets in at least several tens of seconds 
after the explosion (when enough external gas is swept by the shock), we con- 
clude that very early afterglow emission can, in principle, come from radiative 
shocks propagating in relatively strong Wolf-Rayet winds (v4^,> few). 

It is remarkable that under these conditions, all the emission shall come from 
a thin shell of thickness ~ tcooic ~ 1 meter (internal shocks) and ~ 100 km 
(external shocks), that is, from the region of main shock compression. This 
region is already well resolved in 2D PIC simulations. Moreover, one can obtain 
the emitted radiation directly from PIC simulations. 
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